Fluorescence quenching in graphene: a fundamental ruler and evidence for transverse 

plasmons 
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Graphene's fluorescence quenching is studied as a function of distance. Transverse decay channels, 
full retardation and graphene-field coupling to all orders are included, extending previous instan- 
taneous results. For neutral graphene, a virtually exact analytical expression for the fluorescence 
yield is derived, valid for arbitrary distances and only based on the fine structure constant a, the 
fluorescent wavelength A, and distance z. Thus graphene's fluorescence quenching measurements 
provide a fundamental distance ruler. For doped graphene and at appropriate energies, the fluores- 
cence yield at large distances is dominated by transverse plasmons, providing a platform for their 
detection. 
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I. INTRODUCTION 

The optical properties of graphene have attracted im- 
mense attention due to the potential applications in the 
field of photonics and optoelectonics 1 -. Absorption, for 
instance, has the universal value ~ na for light in the vis- 
ible spectrum, depending on the fine structure constant 
a, but not on material's properties^* 3 ., providing a deci- 
sive clue for its original identification-. The large intrin- 
sic carrier mobilities and doping tunability have led to a 
number of proposals for optoelectronic applications 1 -^— , 
where the engineering of long-lived graphene's plasmons 
could play a major role^r— . 

Graphene (and its oxide) exhibits excellent quenching 
of nearby fluorescent materials 1 — ; a property shared 
with carbon nanotube a 15 i 16 . This technique has al- 
lowed spectacular contrast images, enabling far eas- 
ier optical identification 1 ^ (and prospects for device 
manipulation^!) of graphene's flakes. Given the ma- 
ture nature of fluorescent microscopy, particularly in 
the biological sciences, their combination with increas- 
ingly available and versatile graphene nanostructures 
could open interesting research directions 1 ^. Further- 
more, graphene's fluorescence quenching has also been 
proposed as a convenient probe of the intrinsic excita- 
tions of doped graphene such as plasmons 1 ^ 

Here, we theoretically address the effect of graphene 
on the fluorescent material's yield as function of the dis- 
tance z within a unifying formalism. By this, we discuss 
for the first time the transverse decay channels known to 
lead to characteristic features only found in graphene. 19 
The process implies non-radiative (Forste r 20 i 21 ) decay of 
the excited dye, with energy transfer to graphene's ex- 
citations. This mechanism is expected to dominate over 
competing charge-transfer processes except, perhaps, in 
near contact situation o 22 ' 23 . Our motivation stems from 
the seemingly paradoxical coexistence in graphene of very 
strong fluorescence quenching and nominally weak cou- 
pling to the electromagnetic field, as judged from the 
absorption results. 



The distance behavior of fluorescence quenching as 
function of the distance z of a dye molecule with re- 
spect to a metal surface was found long ago to be 
z -4 24 rpkjg i aw repeatedly been found in subse- 
quent studies of energy transfer involving various kinds 
of nanomaterialsj2£ For graphene, the topic has been dis- 
cussed by Swathi and Sebastian, treating the (instanta- 
neous) longitudinal coupling between graphene and fluo- 
rescent material to lowest order, and again obtaining the 
z~ 4 -law. 26 Here, we extend their calculations in three sig- 
nificant ways: i) we include transverse decay channels in 
the calculations ii) the coupling between graphene and 
the electromagnetic field is taken to all orders iii) we use 
the full (retarded) photon propagator. 

Our results for undoped graphene are as follows: i) 
a compact, analytical approximation is obtained and 
shown to provide a virtually exact description of the fluo- 
rescent yield for all distances. This analytical expression 
only depends on the ratio of the distance to the emitting 
wavelength (z/X) and the fine structure constant a. A 
fluorescence measurement thus becomes a distance ruler, 
based only on fundamental constants - a long sought goal 
of the field, ii) Retardation is shown to modify the z -4 
law into a slower z~ 2 behavior, with the transverse decay 
channel dominating at long distances. 

For doped graphene, we show that the fluorescence 
yield is mostly determined by the plasmonic modes 18 
where at certain frequencies transverse plasmons yield 
the dominant, extremely slowly decaying contribution at 
large distances. This should help to detect these modes 
characteristic to graphene and intimately linked to the 
chirality of its elementary excitations. 

The paper is organized as follows. In section II, we 
introduce the formalism defining the atom-field coupling 
in the presence of graphene and the induced fluorescence 
quenching. In section III, we present our results and 
then close with a summary. In an appendix, details on 
the explicit form of the decay rates are given. 
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II. FORMALISM 
A. Atom-field coupling in the presence of graphene 

We consider a fluorescent atom modeled by a two level 
system and the electromagnetic field described in a gauge 
without scalar potential. Within the standard dipole ap- 
proximation, an excited atom will decay into the ground 
state at a rate (7) given by>2£ 



fc 7 = -2uj 2 p* a kiiD a P(r,r;Lj) Pp , 



(1) 



where p is the dipole matrix element between ground and 
excited states, separated in energy by fwj. V represents 
the retarded photon Green's function defined as usual, 



V al3 (r,r' 



dte l ^([A a (r,t),A (r',O)}), 



(2) 

where A(r) is the vector potential operator at the atom's 
location. Notice that Eq. © corresponds to the photon 
propagator without the fluorescent atom, but in the pres- 
ence of graphene. 

In vacuum (that is, without graphene), the pho- 
ton Green's function (r — r';w) is given by the 



Fourier transform of V (k, u>) 



j 2 — c 2 k 2 



(Suf3 



k 2 



2(.qui' 2 ~ S W' • Later inclusion of graphene, assumed per- 
pendicular to the z axis, will preserve the parallel (to 
graphene's plane) component of momentum, q — (gi, 52), 
as a good quantum number. Therefore, it is convenient 
to employ the following representation for the vacuum 
Green's function 

Vf( Z) z';q,oj) = -L J dk z e ik ^-^Vf(k,oj), (3) 

with k — (q,k z ). Physically, Eq. ([3]) represents the 
vector potential in a plane perpendicular to the z-axis 
located at the position z due to currents in a parallel 
plane at location z' . 

The in-plane components of the tensor T>^ , decom- 
posed into longitudinal and transverse contributions, are 
given by 

Vi j (z,z') = d ie -^- z '\ ^ + d t ^- z '\ (Sij - 

q (4) 

with = 1,2 and q' = \J q 2 — (uj/c) 2 . The functions 
di t t(q,uj) are given by (dependencies removed for clarity) 



The presence of a graphene plane at the location z\ 
modifies the vacuum Green's function as follows 

V a \z, z') = Vf(z, z') + 2? m (z, z^x^izuz'), 

(8) 

where sum over repeated indexes is assumed. e is 
the electron charge and X^il: w ) represents graphene's 
current-current total response to external fields. The lat- 
ter, decomposed into longitudinal and transverse contri- 
butions, is given by 



Xi QiQj 



Xi 



1 -<*d lXl q - 1 • x t iSiJ Q 2 ^' (9) 



where we take the non-interacting (RPA) , well-known ex- 
pression for the longitudinal 2 ^ and transverse^ compo- 
nents at zero doping: 



Xi 



9s9v 



16fo yjq 2 ~ {io/v f 



Xt 



9s9v 



v\/q 2 - (uj/v) 2 , 
(10) 

with spin and valley degeneracies, g s = g v = 2, and 
graphene's velocity v. For finite doping, we refer to the 
expressions given in Refsi 29 ' 30 . 

The previous calculation of Swathi and Sebastian 2 ^ 
would correspond to zero doping and retaining only the 
(numerator of the) longitudinal response (x,) in Eq. (0), 
while setting c — > 00 in the photon propagator (instanta- 
neous limit). 



B. Fluorescence quenching 

Consider the graphene sheet placed at the origin (z\ = 
0) and the excited atom at a distance z. The expression 
([1]) for the decay rate can be decomposed as 



h 7nj = ~2u 2 p* a { 



(2tt) S 



d 2 qlmD a!3 (z,z)} pp , 



(11) 

with the (graphene's modified) photon Green's function 
given by Eq. (|8]). The q label classifies the final field 
states into evanescent excitations (q > oj/c), and prop- 
agating excitations (q < oj/c), the latter being the ob- 
served photons. Therefore, the total decay rate is given 
by the radiative and non-radiative contributions to the 
decay rate, 



d t = —. 



Vl*{z,z') 



T S(z 



(5) 



2eow 2 ' ' 2eog' 
The remaining tensor components are written as 

V$(z, z') = V$(z, z') = ^d ; e-«>- 2 'l sgn(z - z') (6) 

and 



(7) 



(12) 



Let us consider the rate of observed photons $. In ad- 
dition to 7rnr , it will depend on the rate at which the 
atom is pumped into the excited state 7 e:rc . Furthermore, 
not all propagating photons are observed, a fraction be- 
ing later absorbed by graphene 7a6s . The excitation rate 
is hardly affected by the presence of graphene and the 
fraction of emitted photons later absorbed is, up to log- 
arithmic corrections, of the order of the fine structure 



3 



constant (a = . e „ ). Therefore, the ratio of the to- 
tal observed fluorescence when the atom is at distance z, 
<&(z), to that at infinite distance, $oo, can be written as 

where the neglected terms amount to minute relative cor- 
rections of order a 2 in the expression (|13p. 

III. RESULTS 
A. Zero doping 

We have evaluated the distance dependence to un- 
doped graphene of the observed fluorescence, Eq. (fT3|). 
with 7 r nr obtained from Eq. (TTT|) and Eq. ((U . There is a 
sharp difference in graphene's effect on the non-radiative 
and radiative contributions to the decay. Graphene mod- 
ification of the (radiative) vacuum decay is a weak effect, 
proportional to a (up to logarithmic corrections). There- 
fore, setting 7 r ~ 7o, with the vacuum decay rate given 
by hr/ = 3 ^. e " c a , the results can be written as 

— « — = /3i7iA + &7 2 / 2 + MiL (14) 



with the physically relevant magnitudes given by 




where the function g(a) can be written in terms of the 
sine and cosine-integrals (si, ci) as 

g(a) = — ci(a) cos(a) — si(a) sin(a). (16) 

The coefficients /3i,2.3 are mere geometric factors de- 
pending on the emitting dipole orientation, with f3\ = 
(p 2 /2+p 2 z )/p 2 , ft = p 2 Jp 2 , and ft = ( P 2 /2)/p 2 . All 
information about graphene in Eq. (ThO) is relegated to 
the dimesionless factors fi, derived and discussed in the 
appendix. 

The first term, 7^, coincides with the unretarded 
contribution to the decay into graphene's longitudinal 
(charged excitations) , previously considered^. The other 
two terms (|15bp and (|15c[) . absent in a non-retarded cal- 
culation, prevail at large distances. The contribution of 
Eq. (|15b[) comes from charged excitations whereas (|15bp 
is due to transverse excitations. Quantitatively, it is the 
last term (|15cp . which provides the dominant large dis- 
tance asymptotic behavior. 
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FIG. 1. Distance dependence of the analytical approxima- 
tion for the decay channels (7, solid lines) compared to 
the exact results (7i/i, dashed lines). Black solid (dashed) 
line: 71 (71/1). Red solid (dashed) line: -y 2 (72/2). Green 
solid (dashed) line: 73 (73/3). (Exact results correspond to 
A = 500nm). Inset: Fluorescence as a function of distance. 
Black line: exact result. Red line: analytical approximation. 

B. Analytical approximation 

Our analytical approximation consists in setting the 
functions fi equal to one. The approximation fi ~ 1 

for i = 1,2 holds when x <C 1 and a;^^ -2 " — 1) 1, 
where xq = y--— , with graphene's fine structure con- 
stant u g = -a. For graphene parameters, the latter 
condition dominates and can be recast as xn < or, 
equivalently, j > j m 10 -3 , justifying our analytical 

approximation. The approximation w 1 applies for 
x < 1, implying f > ^ « 10" 4 . 

The exact and approximate non-radiative decays are 
plotted in Fig. [T] One sees that the approximation 
only fails in the extreme sub-wavelength regime. Fur- 
thermore, even though the exact decay rate saturates for 
(z/X) — > whereas the approximate one diverges, this 
saturation value is so huge that the difference between 
exact and approximate results has virtually no impact 
on $ as seen in the inset of Fig. Q] where the exact and 
approximate curves are indistinguishable. Fig. Q] con- 
firms that retardation channels, 7 2 3 , control the decay 
at large distances, leading to a z~ 2 behavior dominated 
by graphene's transverse excitations of Eq. (|15cl) . 

C. Fundamental ruler 

Notice that only a, z and A appear in expressions (|15l) . 
without any reference to graphene's properties. This im- 
plies that a measurement of the fluorescence quenching 
amounts to a measurement of the distance z, in terms of 
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the light's wavelength A, and the fine structure constant 
a. In others words, it provides us with a fundamental 
distance ruler. 

In general, one would expect graphene's properties to 
drop out from dimensionless optical properties involving 
graphene's excitation within the light-cone, such as the 
absorption. There, u> 3> vq, and graphene's response be- 
comes local, leaving a as the sole coupling scale. But a 
cursory application of this reasoning to our case would 
justify expressions like Eqs. (fT5|) only for distances z > A, 
where graphene's light-cone excitations dominate. Sur- 
prisingly, the analytical expression apply for virtually ar- 
bitrary short distances, implying graphene's excitations 
deep into the evanescent region where the g-dependence 
does not seem obviously negligible. A further surprise 
is the enormous efficacy of graphene as a fluorescence 
quencher, particularly in view of the nominally weak cou- 
pling with the field, set by a. Both facts, range and 
strength, can be explained if the naive nominal range for 
the expected disappearance of graphene's Fermi velocity 
v in the expressions z > A, can be extended to much 
shorter distances z«A. Our approximations show that 
this is indeed the case. The blow up of (A/2:) 4 in the first 
term of Eq. (|14l) when (z/X) — > 0, more than compen- 
sates the overall factor a, leading to strong quenching, 
as experimentally observed 1 - 1 -"— . 



D. Finite doping 

Graphene's fluorescence quenching in doped graphene 
due to graphene's longitudinal coupling to the light field 
has first been analyzed in Reff 1 ^ within the unretarded 
approximation. For large frequencies, Hlj>2Ef, the flu- 
orescence yield of undoped graphene is obtained. For 
frequencies huj^Ep, the quenching behavior is domi- 
nated by longitudinal plasmon excitations, leading to an 
characteristic exponential decay with distance. 

Here, we extend the discussion by also analyzing the 
transverse plasmon excitations. 1 - These exist in the range 
1.667 < 0<2 with SI — uj/(vkp), leading to a decay rate, 
7 t , which dominates at long-distances: 



1 



1 



(17) 



q 2 — (uj/c) 2 , with plasmon momentum 



where q' p = 

q p . Approximating graphene's transverse response by the 
long-wavelength limit, we obtain the following analytical 
expression: 



2l 

7o 



3irp 2 a 



-z/z 



with z, 







4a(vF I c)kF f (fi) and 



In 



2 + n 



2-0 



- 1. 



(18) 



(19) 




FIG. 2. Distance dependence for the decay channels 7^ for 
doped graphene for S7 = 1.75. Black solid line: Red 
solid line: 72. Green solid line: 73. Green dashed line: 7 t . 
Inset: crossover length as a function of frequency. (Numerical 
values correspond to A = 500nm and electronic density n = 



In Fig. [2 the distance dependence of the various decay 
channels is shown for $1 = 1.75 as obtained numerically. 
The response is controlled by the singular features in the 
dispersion relation, leading to an exponential distance 
law. The transverse decay channel 73 is almost entirely 
due to the transverse plasmon mode 7 t , and dominates 
the long-distance behavior beyond a crossover length, z c , 
whose frequency dependence is shown in the inset of Fig. 
[2] We finally note, that the extremely slow decay rate of 
7 t as well as the large distance required for the onset of 
the power law in Eq. (|15c[) can be linked to the condition 
1 - dtxt « 0. 



IV. SUMMARY 

We first studied the fluorescence quenching efficacy 
of undoped graphene as function of distance, including 
transverse decay channels, retardation and graphene-field 
coupling to all orders. For shorter distances, we confirm 
the validity of previous lowest-order, unretarded results, 
albeit with modifications in range and saturation value. 
Retardation changes the distance law to (A/z) 2 , with 
both longitudinal (charged) and transverse (chargeless) 
graphene's excitation contributing to it, the latter dom- 
inating in the truly large distance asymptotic regime. 

A compact, virtually exact analytical expression has 
been obtained for the zero-doping fluorescence yield for 
all distances. It involves only a and z/X, and, therefore, 
endows graphene's fluorescence quenching measurements 
with the unique status of a fundamental ruler. 

We also showed that measurements of the fluorescence 
quenching efficacy of doped graphene at appropriate fre- 
quencies and large distances give direct evidence of the 
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existence of transverse plasmons unique to two-band ma- 
terials like graphene. 

Finally, we note that our results might also be im- 
portant because non-radiative decay of fluorescent ma- 
terials amounts to electron-hole (carrier) generation in 
graphene. Direct (propagating field) photo-generation 
of carriers is inefficient in graphene (absorption w not), 
whereas indirect photo-generation by the evanescent field 
(that is, fluorescence quenching) can be very effective, ul- 
timately controlled by the dye's absorption and quantum 
yield. This suggests new ways of enhancing graphene's 
photo-responsivity^Ir— . 
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where xo = j--—, with graphene's fine structure con- 
stant a g — ^a. The approximation fx tzs 1 holds when 

io < 1 and £o(~~T i — -0 ^ 1- For graphene parame- 
ters, the latter condition dominates and can be recast as 
x < or, equivalently, f > § ~ 10~ 3 , justifying our 

approximation /i ~ 1. 

Now, we consider the non-radiative decay channel cou- 
pling graphene's longitudinal response with an out-of- 
plane dipole. Using Eqs. (JSJ) and ©, Eq. (|A1|) leads 
to 



c 2 q' 2 



Im- 



dfxt 



1 - d iXi 
(A5) 

The integration can be rewritten as the two first terms 
of Inrllr given in the main text, where the new dimen- 
sionless factor fi, corresponding to the second term in 
the sum of Eq. (|A5I) . appears. It is given by 



Appendix A: Explicit form of decay rates 

Here we outline the derivation of the non-radiative de- 
cay channels at zero doping 



2J 2 p* a { 7 ^ [ d 2 qImD^(z,z)}p . 

( 27T ) Jq>u/c 



(Al) 

As the non-radiative decay rate only exists in the pres- 
ence of graphene, we classify its contributions accord- 
ing to the nature of graphene's excitations: longitudi- 
nal (charged) and transverse. Longitudinal (I) excita- 
tions couple to both in-plane (p M ) and out-of-plane (p z ) 
components of the dipole matrix elements. Graphene's 
transverse (t) excitations only couple to the in-plane (p u ) 
dipole matrix element. Therefore, we write for the non- 
radiative decay channels 

7„r = 7j,„ + %z + 7t,„ • (A2) 
Let us consider j t M first. Selecting the longitudinal 
part of Eq. (|9]) in Eq. ([8j , it is straightforward to show 
that Eq. ([AT]) leads to 



-2ty 



POO 

/ q'dq'e~ 2q ' z Im- 



(A3) 



4tt 7 * 1 -dix h 

where we have traded q for q' in the integration. Notice 
that, although the exact (i.e. retarded) photon propaga- 
tors are used, the only manifestation of light's velocity c 
in 7 ; M would be the substitution: v~ 2 — > v~ 2 = v~ 2 —c~ 2 , 
within the square root of Xi at zer0 doping. The quanti- 
tative irrelevance of this replacement makes the instan- 
taneous approximation for r ) l n correct. The integration 
in Eq. (|A3|) then leads to the first term of ^nr/lr given 
in the main text, with the dimensionless factor f\ given 
by 

1 f 1 x 3 e~ x / x °\ / l — x 2 
■h = — I dx 1 ^7^—, (A4) 



3\x£ 



L o Jo 



dx ■ 



(A6) 



The approximation / 2 ~ 1 has the same range of validity 
as that of f\. 

Finally, we consider graphene's transverse excitation 
channels. Using the transverse components of Eqs. Q 
and © in ©, Eq. (fM]) leads to 



h lt u = -2w 2 |^- J q'dq'c- 2q ' z Imj 



d 2 t Xt 
- dtXt ' 



(A7) 



an integration that can be rewritten as third term of 
lnr/lr given in the main text, with the corresponding 
dimensionless parameter /3 given by 



h 



i 



dx ■ 



l(27r 2 a x )J x i + (%a) 2 (l-x 2 y 
with the function g(a) defined^ as 



(A8) 



5(a) 



dx 



xe 



-ci(a) cos(a) — si(a) sin(a). 

(A9) 



The approximation /3 ~ 1 applies for xq <C 1, implying 
v > -r- ~ 10 -4 , as stated in the main text. 

A ~ 47TC 5 

The function g(a) exhibits the following asymptotical 



1 „ 3! 5! 
g(a»l) = — (1 - — + — -...) 



(A10) 



g(a < 1) = -log(e L a) 



with Euler constant C. The first limit in Eq. (|A10I) for 
a = 2ir 2 aj 1 explains the large distance asymptotic 
behavior of j 3 as discussed in the main text. 
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